New solutions of charged regular black holes and their stability 
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We construct new regular black hole solutions by matching the de Sitter solution and the Reissner- 
Nordstrom solution with a timelike thin shell. The thin shell is assumed to have mass but no pressure 
and obeys an equation of motion derived from Israel's junction conditions. By investigating the 
equation of motion for the shell, we obtain stationary solutions of charged regular black holes and 
examine stability of the solutions. Stationary solutions are found in limited ranges of 0.87L < m < 
1.99L, and they are stable against small radial displacement of the shell with fixed values of m,M, 
and Q if M > 0, where L is the de Sitter horizon radius, m the black hole mass, M the proper 
mass of the shell and Q the black hole charge. All the solutions obtained are highly charged in the 
sense of Q/m > 2\/3 £s 0.866. By taking the massless limit of the shell in the present regular black 
hole solutions, we obtain the charged regular black hole with a massless shell obtained by Lemos 
and Zanchin and investigate stability of the solutions. It is found that Lemos and Zanchin's regular 
black hole solutions given by the massless limit of the present regular black hole solutions permit 
stable solutions, which are obtained by the limit of M — > 0. 

PACS numbers: 04.70.-s 



I. INTRODUCTION 

One of the most interesting questions in general rela- 
tivity concerns the inner structure of black holes. How- 
ever, it is hard to give a definite answer to this question 
because of the existence of a spacetime singularity where 
the curvature diverges indefinitely and general relativ- 
ity breaks down. Penrose and Hawking proved that the 
gravitational collapse with physically reasonable initial 
conditions inevitably leads to the formation of a singu- 
larity, known as the singularity theorems [H-Q ■ Although 
such singularities are supposed to be concealed by the 
event horizon and to be isolated from the domain of pre- 
dictability, this means that we cannot describe the entire 
spacetime by the present physics. However, there may be 
situations where some of the assumptions in the singu- 
larity theorem may not be applied, such as the existence 
of the cosmological constant somewhere in the spacetime 
region. Thus, it is of interest if we can construct some 
models of black holes without spacetime singularities. 

Black holes having the regular centers are called reg- 
ular black holes or nonsingular black holes. Existing 
regular black hole solutions may be divided into two 
classes. Solutions belonging to one class are character- 
ized by the property that the black hole spacetime is suf- 
ficiently smooth everywhere. Bardeen gave this type of a 
solution for the first time |4| . The metric of his solution 
asymptotically approaches the de Sitter and the Reissner- 
Nordstrom solutions in the limit of r — > and r — >• oo, 
respectively, where r is a Schwarzschild-type radial coor- 
dinate. If one chooses appropriate parameters for the 
Bardeen solution, its spacetime has two horizons and 



looks like a Reissner- Nordstrom black hole with a regular 
center. Although it was thought that the Bardeen solu- 
tion cannot be an exact solution of Einstein equations, 
Ayon-Beato and Garcia [9] showed that the Bardeen solu- 
tion is given as a gravitational field coupled to a nonlinear 
magnetic monopole. So far, there have been many inves- 
tigations based on Bardeen's work for uncharged cases 
0-0| and for charged cases f§-fi~lj|. The other class of 
regular black hole solutions is composed of the solutions 
constructed by matching two distinct spacetimes with 
a thin transition layer or surface. Typical solutions of 
this class are composed of a single regular de Sitter core 
and exterior black hole spacetime between which a sin- 
gle thin shell exists |13 - [l4j . The layer, which must be 
located within the event horizon because we consider the 
regular black hole, can be either a spacelike dHHH, time- 
like [3] or null hypersurface [l4[ • The regular uncharged 
spherically symmetric black holes [l2j are motivated by 
the assumption that the spacetime curvature has an up- 
per limit which is of the order of the Planck scale and 
the quantum effects become dominant so that the forma- 
tion of the singularity is avoided. The collapsing matter 
will turn into a de Sitter phase when the curvature ap- 
proaches a critical value. This idea was first suggested 
for the cosmological context by Sakharov [16J and Glincr 

As mentioned before, studies on the regular black holes 
are closely related to matching problems of two different 
spacetimes and motion of the thin shell. Israel derived 
convenient matching conditions of spacetimes for non- 
null transition layers [l8j]. Barrabes and Israel general- 
ized Israel's junction conditions to a unified description 
including null hypersurface cases 19 1. (For the spher- 
ically symmetric cases, see Ref. 20|.) As for motion 
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of the shell, gravitational collapse of a charged shell has 
been studied in Refs. [2ll424| . (For the cases of higher 
dimensional spacetimes and other gravitational theory, 
see, e.g., Refs. [IflHil- ) I R those studies, the interior 
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spacetime of the shell is assumed to be fiat [22J, |24f or 
the Reissner-Nordstrom solution with different mass and 
charge from outside one As argued in those studies, 
the shell can be stationary and stable only if the shell has 
pressure. In cases of no pressure, the shell keeps on col- 
lapsing or expanding, or the shell collapsing (expanding) 
at the beginning will turn to expand (collapse). 

Stability of regular black holes is also important be- 
cause unstable solutions cannot occur in nature. For the 
regular black holes with shells, instability of the station- 
ary shell immediately implies instability of the regular 
black hole. Balbinot and Poisson jl5j have analyzed sta- 
bility of spacelike shells of the uncharged spherically sym- 
metric regular black holes considered by Frolov et al. [12j . 
They showed that for a certain parameter, the shell can 
be stationary and stable. In this study, we will apply 
Balbinot and Poisson's method to the cases of charged 
regular black holes; the regular black holes consist of a 
timelike massive charged shell which separates the de Sit- 
ter and the Reissner-Nordstrom spacetime inside the in- 
ner horizon of the Reissner-Nordstrom solution. Lemos 
and Zanchin have considered this type of charged regu- 
lar black holes and obtained exact solutions assuming the 
shell is massless and pressureless [l4j . 

Our aim in this study is threefold; to find new regular 
black hole solutions having a massive thin shell, to ex- 
amine their stability, and to examine stability of Lemos 
and Zanchin's regular black holes. For simplicity, we as- 
sume the shell is constructed of dust, i.e., the shell has 
mass but no pressure. Although the shell is pressureless 
in this study, we consider the de Sitter spacetime inside 
the shell, i.e., there exists matter that corresponds to a 
cosmological constant inside the shell. Thus, the pres- 
sureless shell can be in stationary states. Balbinot and 
Poisson [r| considered that the de Sitter horizon is of the 
order of the Planck scale and it is much smaller than the 
event horizon. However, we do not a priori make any as- 
sumptions for physical scales of the parameters regarding 
the regular black hole in this study. 

The plan of this paper is the following. In Sec. II, we 
briefly describe formalism for a thin shell using the 3+1 
decomposition of Einstein equations and derive equations 
of motion for a thin dust shell. In Sec. Ill, we show 
results of new regular black hole solutions and their sta- 
bility. Stability of Lemos and Zanchin's regular black 
holes is argued in Sec. IV. Then, the conclusion is in the 
last section. Throughout this paper, we use the units of 
c = G = 1, where c and G are the speed of light and the 
gravitational constant, respectively. 



II. FORMULATION 

A. Preliminary 

As mentioned before, we consider solutions of Einstein 
equations in which two different exact solutions, the de 
Sitter and Reissner-Nordstrom solutions, are matching 



by a massive thin shell. Following Ref . [l5| . in this sub- 
section, we concisely describe the formalism treating mo- 
tion of the thin shell sandwiched between two arbitrary 
solutions. 

Let V be the four dimensional spacetimes that have 
metrics g^g and a system of coordinates x± a . Let 
S be a hypersurface described by intrinsic coordinates 
£ a = (C 1 ^ 2 ^ 3 ) an d located at the boundaries of V + 
and V~ . Here and henceforth, we use the greek and the 
roman lowercase letters to describe indices of the four- 
dimensional spacetime and of the three-dimensional hy- 
persurface, respectively. Let n a be a unit normal vector 
to the hypersurface S. Thus, n a has to satisfy 

n a n a = e, e> Q = 0, (2.1) 

where e" is the basis vector on S, defined by 

„ dx a 



(2.2) 



Here, e = 1 (e = — 1) when the hypersurface is timclike 
(spacelike). The induced metric h ab and the extrinsic 
curvature K ab associated with E are, respectively, defined 
by 



h a b = g a pe"e b , K ab 



a /3 



(2.3) 



Here and henceforth, we denote the covariant differen- 
tiation associated with g a p and h a b by the stroke (|) 
and the semicolon (;), respectively. To describe mo- 
tion of the three-dimensional hypersurface, it is useful 
to rewrite the basic equations in the three-dimensional 
form. These equations can be derived by contracting the 
four-dimensional quantities by n a and/or e". By using 
the Einstein tensor G a p contracted by n a and/or e", and 
the Gauss-Codazzi equations, we obtain 



- 2eG afs n a n fi = 3 R + e(K ab K ab 
G a pe a y = K, a - K b a . tb , 



K 2 



(2.4) 



where 3 R is the three-dimensional Ricci scalar associated 
with h ab and K = h ab K ab . The energy- momentum ten- 
sor on the hypersurface, S ab , is given by a jump of the 
extrinsic curvature on £ (see, e.g., Ref. jig). 



tTrS ab = e([K ab ] - h ab [K}), 



(2.5) 



where [K ab ] — K ab \ + — K ab \~ . is evaluated on £ 

by taking limit from V . Then, the energy-momentum 
conservation equation on the hypersurface will yield 



SI 



b:a 



(2.6) 



So far, we have shown the energy-momentum conser- 
vation equation in the general form. We next show how 
these equations are given for the dust shell. If the shell 
is composed of dust, the stress-energy tensor of the shell 
is given by 



S a b 



au a u bl 



(2.7) 
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where a is the surface energy density of the shell and 
u a is the matter velocity on S if the shell is a timelike 
hypersurface. Thus, the energy-momentum conservation 
leads to 



K):a= [T aP U a n fi ], 



(2.8) 



where u a is the four- velocity of the shell given by u a = 
u a e" . The transverse acceleration of the shell is expressed 
by 



a a = u^u 13 = u a h u h e a a + eu a u b K ab n c 



(2.9) 



We are interested in the normal component of a a , which 
describes the motion of the shell. It is straightforward to 
see that [n Q a a ] = u a u b [K ab ]. Equation (12.51) is equivalent 
to 

[K ab ] = 87re (s ab - |/i a6 Y with S = h ab S ab . 

(2-10) 

Then, we have an equation of motion of the shell, 

n a a a \ + - n a a a \- = inca . (2.11) 

We may also construct an equation from the arithmetic 
means of the extrinsic curvatures, (K* b + K~ b )/2. How- 
ever, it is not a useful equation for the present situation. 

B. Equation of motion of the shell 

In order to have a regular center by matching the 
de Sitter and the black hole spacetimes with a mas- 
sive thin shell, at least two horizons including extremal 
cases, in which two horizons coincide, are required [27| . 
For uncharged spherically symmetric cases, since the 
Schwarzschild black hole, which is the outer solution for 
this situation, has a single horizon (event horizon), the 
second horizon must be the de Sitter one. This implies 
that the shell has to be located between the outer event 
horizon and the inner de Sitter horizon and that it nec- 
essarily has to be spacelike. For charged and/or rotat- 
ing cases, the black hole solution has double horizons, 
the event and Cauchy (inner) horizons. We may, there- 
fore, choose any type of shell — timelike, spacelike, or null. 
Since it is physically natural to assume that the shell is a 
timelike hypersurface, we choose the shell to be located 
inside the inner horizon of the black hole solution in these 
cases. 

Assuming V + and V~ to be the Reissner-Nordstrom 
and de Sitter spacetimes, respectively, we apply the for- 
malism given in the previous subsection to the present 
situation. We derive the equation of motion for the shell, 
given by r = R(t), from Eqs. (j2~g)) and (|2~TT1) . with r 
being the proper time of the shell. In this study, as men- 
tioned, we assume the shell to be a timelike hypersurface. 
Thus, the shell radius R is assumed to satisfy R(t) < r_ 
and R(t) < L, where r_ and L denote radii of the inner 
horizon of the Reissner-Nordstrom solution and the de 



Sitter horizon, respectively. The spherically symmetric 
metric that expresses the inside and outside of the shell 
is written by 

ds 2 = - f (r)dt 2 + — !— dr 2 + r 2 (d9 2 + sin 2 9 d(/> 2 ), (2.12) 
and the function f(r) is given by 



f(r) 



f dS {r) = 1 - — , for r < R(t), 
1m (l 2 

f RN (r)^l- — + ^, for r>R(r), 



(2.13) 

where m and Q are the mass and the charge of the 
black hole, respectively. Here and henceforth, Hds 
and Hrn mean functions H evaluated by the de Sit- 
ter and the Reissner-Nordstrom solutions on S, respec- 
tively. Due to the assumptions, R(t) < r_ and R(t) < L, 
fds( r ) is always positive and /ijAr(r) has two roots (one 
root) for Q < m (Q = m). The electric potential is 



Q 



dt and A^dx-* 1 = 0. The four veloc- 



ity of the shell is u a = dx a /dr = (i, R, 0,0), where the 
dot denotes the derivation with respect to r. Because of 
u a u a = — 1, t may be written as 



R 2 + fir) _ p 



f(r) ~ f(r) ' 



(2.14) 



where we have set that t increases with t. The normal 
vector to the shell is, from u a n a = and n a n a = 1, 
given by 



R 

W) 



,13,0,0), 



(2.15) 



where we have considered the normal vector pointing 
from the de Sitter spacetime to the Reissner-Nordstrom 
spacetime. The induced metric on the shell is 

(ds 2 h = -dr 2 + R 2 (T)(d6 2 + sin 2 6 d<j> 2 ). (2.16) 

Since the four-velocity u a and the normal vector n a do 
not have 9- and (^-components, we only need to consider 
t and r components of the T Q ^, which leads to 



dS 



3ffq,3 
87TL 2 ' 



(T aP ) RN = -^(d r A t ) 2 , for a,/3 = 0,1. 

07T 



(2.17) 



Then, the energy conservation equation, 

(au a ), a = [T aP u a n fi ] =0, (2.18) 

leads to 

(i? 2 cr)- = 0. (2.19) 
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Note that T a pu a n" vanishes when T a p oc g a p. If we de- 
fine the proper mass of the shell by M = 4irR 2 a, the 
above equation means M is independent of r. Nonva- 
nishing components of the extrinsic curvature are given 
by 

n a a a =K TT = i, K% = K+ = -L (2.20) 
R it 

From Eq. (12+11) , thus, we have 

pRN~PiS = ^Ra = -4n(Ray, (2.21) 

where we have used (R 2 <r)' = 0. Integrating Eg. (12. 211) . 
we obtain 

r~ w nz 2^ M 

v R2+i -v-v R2+i -ir + w=R +c ' 

(2.22) 

where C is an integration constant. From Eqs. (|2.5[) and 
we have, -{P RN -Pds)/R = [K 6 e ] = [K*,] = 4tt<j. 
Thus, the integration constant C appearing in Eq. (I2.22[) 
has to be zero. Rather than using Eq. ([2.22j ). we employ 
a more convenient form, 

R 2 + V(R) = -1, (2.23) 

where 

V{R) = _(£ + %-^_K)\*. (2.24) 

17 I 2M 2RJ L 2 K ' 

Equation (|2.23p is a kind of energy conservation law be- 
cause one may interpret V(R) as an effective potential. 
This equation is nothing but an equation of motion for 
the massive thin shell located on the surface of the de 
Sitter sphere. The stationary solutions, R = const., can 
be obtained by solving V(R) = -1 and dV{R)/dR = 
simultaneously. The stability of the stationary solutions 
can be checked by the condition d 2 V (R) / dR 2 > at the 
stationary point, i.e., whether or not the stationary shell 
is at a local minimum of the effective potential V(R). 
From Eq. (j2~2"4"f . we see V = V(\M\). To obtain values 
of M, we use Eq. (|2.22j) after obtaining solutions of the 
regular black hole. Since, in this study, we are concerned 
with the regular black hole solution, we assume that no 
naked singularity occurs, i.e., m > Q. 

III. RESULTS 

In order to show numerical results, we employ the units 
of L = 1, e.g., R/L -> R, M/L -> M, m/L -> m, and 
Q/L — » Q. Let us describe a method to obtain the equi- 
librium states of the regular black hole. To obtain sta- 
tionary solutions numerically, we solve V(R) + 1 = and 
dV (R) I dR = simultaneously with a Newton- Raphson- 
like iterative scheme. During the iteration procedure, 
values of \M\ and m are kept constant. Thus, the two 



algebraic equations, V(R) + 1 = and dV(R)/dR = 0, 
contain only two unknown parameters, R and Q. Then, 
we can obtain a regular black hole solution if the iter- 
ation procedure successfully ends. After obtaining so- 
lutions, we obtain values of M from Eq. (|2.22l) and 
check their sign of d 2 V(R)/dR 2 to see their stability. If 
d 2 V{R)/dR 2 > {d 2 V(R)/dR 2 < 0), then the solution is 
stable (unstable). Since, as mentioned before, we assume 
that there is no naked singularity, and the shell is inside 
the inner horizon of the black hole solution, we are con- 
cerned with the solutions satisfying conditions m > Q, 
R < 1, and R < r_. Otherwise, we do not admit the so- 
lutions as those of the regular black hole model that we 
consider in this study. Although the solutions of M < 
are not physically acceptable in normal situations, they 
are allowed in Eq. (|2.23[) and might be useful in some ex- 
otic situations. Thus, we show the results of the M < 
case as well in this study (see, e.g., Ref.[24j]). 

It is helpful to examine properties of the effective po- 
tential V in order to understand how the regular black 
hole solution is obtained. In Figs. 1 and 2 and Figs. 3 
and 4, we show typical behaviors of the potential as func- 
tions of R for stable and unstable stationary solutions, 
respectively. Figures 2 and 4 are magnified figures of 
Figs. 1 and 3 around extremal points, respectively. The 
potential of the stable (unstable) configuration, given 
in Figs. 1 and 2 (Figs. 3 and 4), are characterized 
by a set of parameters (m,Q,\M\) = (1.2,1.191,0.132) 
[(m,Q,\M\) = (1.31,1.309,0.01)]. A local minimum 
(maximum) of the potential, shown in Figs. 1 and 2 
(Figs. 3 and 4), is at R = .87661 (R = 0.86623). As 
can be seen from Eq. (I2.24[) , the effective potential V 
diverges to minus infinity as R — > and as R — > oo, 
which means that V has at least one maximum. This 
property may partly be confirmed in Figs. 1-4. It is also 
observed in Figs. 2 and 4 that V « — 1 at the extremal 
points. Thus, one sees that these potentials permit the 
stationary solutions. 

To investigate basic properties of the regular black hole 
solutions, we calculate many sequences of stationary solu- 
tions. The sequences of stationary solutions are specified 
by a fixed parameter m, and are obtained by increasing 
a parameter \M\ from a minimum value 10 -4 (see, also, 
the first paragraph of this section). In other words, they 
are given as a set of functions R(M) and Q(M) with a 
fixed parameter m. The sequences of stationary regular 
black hole solutions, given as functions R(M) and Q(M), 
are shown in Figs. 5 and 6. As shown in these figures, we 
obtain the regular black hole solutions for positive and 
negative values of M. One of the interesting findings in 
this study is that all the positive (negative) M solutions 
are stable (unstable), and it seems that these stable and 
unstable solutions belong to continuous sequences of the 
stationary solutions of a fixed m as can be seen in Figs. 
5 and 6. For the sequences of the solutions obtained in 
this study, R's are increasing functions of M (see Fig. 

5) , and (Q/m)'s are decreasing functions of M (see Fig. 

6) . For high m sequences of the solutions, values of R 
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approach unity (see Fig. 5). It is important to point 
out that all the solutions obtained in this study satisfy 
1 > Q/m > 0.86, i.e., they are highly charged, and that 
some solutions are nearly extreme in the sense Q ~ to at 
one end of the sequences of the solution (see Fig. 6). 

As shown in our numerical results argued so far, the 
regular black hole solutions are only found in some re- 
stricted parameter regions. Thus, it is useful to show 
clearly in which parameter regions the regular black hole 
solutions with a timelike thin shell occur. In Figs. 7 and 
8, we show two-dimensional parameter regions where reg- 
ular black hole solutions are found in this study. Stable 
and unstable solutions are found in 0.87 < m < 1.99 with 
positive M and in 1.3 < m with negative M, respectively 
(see Fig. 7). As can be observed in Fig. 7, the maxi- 
mum value of M (~ 0.43) for the solutions obtained in 
this study is achieved by the solution with the minimum 
value of to (« 0.87). From Fig. 8, it is observed that 
the radius of the shell is in the range 0.866 < R < 1 for 
the solutions obtained. In the present study, we find no 
upper limit of to for existing unstable regular black hole 
solutions. Note that we show results of unstable solu- 
tions only for the range of 1.3 < to < 1.99 in Fig. 7, but 
to = 1.99 does not mean the upper limit. For the cases of 
to > 2, we cannot determine the minimum value of \M\ 
for unstable solutions because of numerical difficulties. 



IV. DISCUSSIONS 

As mentioned in the Introduction, one of our aims in 
this study is to analyze the stability of the charged reg- 
ular black holes with a massless thin shell derived by 
Lemos and Zanchin |Tij . In order to investigate regular 
black holes with a massless thin shell, we may consider 
the limit of M — > in our formalism. We focus on the 
M > case in this section since the M < case is not 
physically acceptable in usual situations. Besides, all the 
solutions with M < are unstable, so they are not fea- 
sible as the regular black hole models. Thus, we may 
exclude the M < case for our physical interests. In 
the M — > limit, the effective potential V and its first 
derivative dV/ dR are approximated by 



V 

dV 
dR 



2 2mR+% 



4M 2 R 2 
Q 2 -3fi) (Q 2 -2TOi?+£± 
2M 2 R 3 



(4.1) 



(4.2) 



Thus, we may obtain stationary solutions in the limit of 
i\/ — ^ if the following conditions are satisfied: 



Q 2 - 2mR O(M) , 



i? 4 



(4.3) 
(4.4) 



We then assume that charged regular black hole solutions 
in the massless limit may be expanded in terms of M as 
follows, 



R 2 = QL/V3 + -ALM + 0(M 2 ) , 
R 3 

m = 2—r + BM + 0{M 2 ), 



(4.5) 
(4.6) 



where A and B are functions independent of M . Substi- 
tuting Eqs. (|4.5p and (|4.6p into V and dV/dR, we obtain 



V 

dV 
dR 



R 2 
L 2 



(2AR + 2BL) 2 



4L 2 

2ABL -2R+ 2A 2 R 
L 2 



+ O(M) , (4.7) 
0{M) . (4.8) 



Then, stationary solutions in the limit of M — > are, in 
terms of Rq and L, given by 



Qo 

A 
B 



V3 2 
— Ro , 

L 2 ' 

Ro 

7^ 

l- A 2 
AL 



R l 



-Ro 



(4.9) 
(4.10) 
(4.11) 

(4.12) 



where quantities indicated by the subscript correspond 
to stationary solutions in the limit of M — > 0. From 
the conditions for the regular black hole with a timelike 
thin shell, r_ > Rq, L > Rq and mo > Qo, we obtain a 
constraint for Rq, given by L > Rq > ^f-L. These mass- 
less limit solutions, given by Eqs. (|4.9|) and (|4.10|) . are 
exactly the same as those given by Lemos and Zanchin, 
although our notation is different form theirs. [Lemos 
and Zanchin also derived the relations R > Q/V3 and 
to < 2Q/v / 3- (These inequalities are derived from Eqs. 
(gU) and (|4TTU1) and L > Rq > For these mass- 

less limit solutions, the second derivative of the effective 
potential d 2 V/dR 2 is approximated by 



d 2 V r URq^/L 2 



R 2 



dR 2 



ML 3 



(4.13) 



where the conditions for the stationary solution (|4.5|) . 
(UU), and (|4.11IH(|4.10p have been used. For the sta- 
tionary solution, we therefore see that d 2 V/dR 2 — > oo as 
M — > 0. This means that in the massless limit M — > 0, 
we have stable solutions. Note that as can be seen from 
Eqs. (|4.1ip and (|4.12p . the massless limit solutions break 
down when Rq — > L. This is because as Rq — > L, the shell 
becomes lightlike, for which the present formalism for the 
timelike shell is not applicable. 

The above results show that in the limit of M — > 0, the 
solutions exist only for v / 3/2(w 0.866025) < R/L < 1, 
3 V / 3/4(« 1.29904) < m/L < 2, and y/3/2 < Q/m < 1. 
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Our numerical solutions with M — 10~ 4 will satisfy the 
massless condition (M = 0) with good accuracy. Thus, 
we may regard these solutions as good approximations for 
exactly massless solutions. Some of those solutions are 
displayed in Table I. In this table, we may confirm that 
infinitesimal quantities R 2 — QLj\fi and m—2R 3 /L 2 and 
a divergent quantity d 2 V/dR 2 indeed depend on M as 
given by Eqs. (|4"3)) . P~o) . and P~T3]) . We see that sets of 
parameter (m, M, R, Q) = (1.3, 10" 4 , 0.86626, 1.29964) 
and (to, M, R, Q) = (1.99, 1CT 4 , 0.99865, 1.72626) corre- 
spond to the lower and the upper endpoints, respectively, 
in the relation of the Lemos and Zanchin solutions (Q <) 
to < 2Q/V3. 

These analyses show that Lemos and Zanchin's regu- 
lar black hole solutions, given by the massless limit of 
the present regular black hole solutions, permit stable 
solutions. This conclusion, however, is not a final one 
because a properly massless thin shell case is excluded in 
the present analysis and because we only consider an ex- 
ample of regular black hole solutions that coincide with 
Lemos and Zanchin's solution in a massless thin shell 
limit. Thus, further analyses are required to draw a 
definite answer of whether Lemos and Zanchin's regular 
black hole solutions are stable or not. 

Comparing these analytic results discussed so far to 
numerical results given in the last section, we may guess 
the lower and upper limits of the physical quantities for 
the regular black hole model to exist for the case of M =^ 
0. Then, it may be conjectured that m/L < 2, \/3/2 < 
Q/m < 1 and V3/2 < R/L < 1 for the stable regular 
black hole solutions (see Figs. 7 and 8). 

Finally, let us consider a physical scale of the stable 
regular black hole solutions we obtain in this study, which 
have not been specified so far. For the stability analysis 
of Schwarzschild-type regular black holes in Refs. [l2l. Il5jj . 
the de Sitter horizon radius is assumed to be of the order 
of the Planck scale and L R < r + , where r + denotes 
the event horizon radius. Thus, L « l p , where l p denotes 
the Planck length. On the other hand, if we take the 
above assumption, i.e., L w l p , our analysis, based on a 
classical approach, fails since the present solutions satisfy 
R < l p ~ L < r_|_. Since the curvature invariant of the de 



Sitter spacetime is R^ V paR tlvf " J = 24/L 4 , if there exists 
an upper bound of the curvature and our analysis is valid, 
the upper bound of the curvature has to be smaller than 
that of the Planck scale. However, our results show that 
if we assume the de Sitter horizon radius is of the order of 
the Planck scale, the present stable charged regular black 
hole solution is restricted to quantum size. Although the 
de Sitter horizon radius is assumed to be other physical 
scales of the vacuum phase transition, such as the grand 
unified theory (GUT) scale, the present stable black holes 
are also restricted to quantum size. 

V. CONCLUSION 

We have constructed new regular black hole solutions 
by matching the de Sitter solution and the Reissner- 
Nordstrom solution with a timelike thin shell. The 
thin shell is assumed to have mass but no pressure and 
obeys an equation of motion derived from Israel's junc- 
tion conditions. By investigating this equation of motion 
for the shell, we obtain stationary solutions of charged 
regular black holes and examine stability of the solu- 
tions. Stationary solutions are found in limited ranges 
of 0.87L < in < 1.99L, and they are stable against small 
radial displacement of the shell with fixed values of to, Q, 
and M if M > 0. All the solutions obtained are highly 
charged in the sense of Q/m > 2\/3 « 0.866. By taking 
the massless limit of the shell in the present regular black 
hole solutions, we obtain the charged regular black hole 
with a massless shell obtained by Lemos and Zanchin [l4| 
and investigate stability of the solutions. It is found that 
Lemos and Zanchin's regular black hole solutions permit 
stable solutions. 
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TABLE I: Physical quantities for the stationary solutions with 

|M| = icr 4 . 



R 


Q 


m 


R 2 - QL/V3 


m - 2R 3 /L 2 


d 2 V/dR 2 


0.866 


1.29 


1.30 


5.78 x 10~ 5 


-1.00 x 10~ 4 


5.19 x 10 4 


0.908 


1.43 


1.50 


7.31 x 10~ 5 


-1.55 x 10~ 4 


4.55 x 10 4 


0.947 


1.55 


1.70 


9.86 x 10~ 5 


-2.48 x 10 -4 


3.63 x 10 4 


0.998 


1.72 


1.99 


6.41 x 10~ 4 


-5.10 x 10~ 4 


5.47 x 10 3 



potential 




0.2 0.4 0.6 0.8 1 



FIG. 1: Effective potential V(R) for the stable solution with 
m = 1.2, Q = 1.1913 and |M| = 0.132. The local minimum 
is at R = 0.87661. 
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potential 
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0.8 0.84 0.88 0.92 0.96 1 

R/L 

FIG. 2: Magnified figure of Fig. 1 around the local minimum. 

potential 




0.6 0.7 0.8 0.9 1 
R/L 

FIG. 3: Effective potential V(R) for the unstable solution 
with m = 1.31, Q = 1.3909 and \M\ = 0.01. The local 
maximum is at R — 0.86623. 
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FIG. 4: Magnified figure of Fig. 3 around the local maximum. 
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=0.9 
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FIG. 5: Radius of the shell R for the sequences of the sta- 
tionary regular black holes with fixed m, given as functions of 
the mass of the shell M. From bottom right to top left, the 
curves correspond to the sequences with m = 0.9, 1.1, 1.3, 
1.5, 1.7 and 1.9. 



relation of M, Q and m 





FIG. 6: Ratio of the black hole charge to the black hole mass 
Q/m for the sequences of the stationary regular black hole 
with fixed m, given as functions of the mass of the shell M. 
From top right to bottom left, the curves correspond to the 
sequences with m = 0.9, 1.1, 1.3, 1.5, 1.7 and 1.9. 
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relation of M and m 




0.8 1 1.2 1.4 1.6 1.8 



FIG. 7: Parameter region in (m, M) plane where the regular 
black hole solutions are found. The stable (unstable) solutions 
exist in the positive (negative) M region as indicated by the 
labeles. 
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0.92 0.94 0.96 0.98 
R 



FIG. 8: Parameter region in (R, Q) plane where the regular 
black hole solutions are found. The stable (unstable) solutions 
exist in the domain bounded the dotted (short dashed) curves. 
The long dashed curve indicates the boundary between the 
stable and unstable solution regions, given as the \M\ = 10~ 4 
sequence of the solutions. From bottom to top, the solid 
curves correspond to the sequences with m = 0.9, 1.1, 1.3, 
1.5, 1.7 and m = 1.9. 



